Theoretical Study of Magnetic Moments Induced by Defects at the SiC (110) Surface 
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The effect of different surface defects on the atomic and electronic structures of cubic /3-SiC(110) 
surface are studied by means of a first principles calculation based on Density Functional Theory 
using the SIESTA code. In the calculations, the possibility of different spin population at each atom 
is allowed. We find that while adsorption of oxygen or nitrogen and adsorption of hydrogen at the 
C surface carbon atoms do not induce magnetic moments on SiC(llO); Si vacancies, substitutional 
C at the Si site and H or F adsorbed at the silicon surface sites induce localized magnetic moments 
as large as 0.7 /j,b at the carbon atoms close to the defect. The local magnetic moment arrangement 
varies from ferromagnetic in the case of H adsorption to antiferromagnetic in the Si vacancy and 
substitutional C cases. The case of H adsorption at the Si surface atoms is discussed in detail. It is 
concluded that magnetism is mainly due to the local character of the C valence orbitals. 



PACS numbers: 71.20.Mq,75.70.Rf, 81.05.Uw 

The localized character of the carbon valence charge 
has motivated the study of the possibility of carbon 
atoms developing localized magnetic moments near de- 
fects PQ. This has prompted the study of a variety 
of defects at different carbon atomic configurations like 
graphite [5] , graphene [31 0] , nanotubes [51 [5] and other 
configurations [71 [H] including defect-free diamond sur- 
faces [9]. Here, in this work, we consider the appear- 
ance of localized magnetic moments near defects on the 
(110) surface of cubic silicon carbide (/?-SiC). The quasi 
one dimensional character of the Si-C chains at the sur- 
face and the localized character of the carbon atomic 
valence charge suggest the possibility of defect induced 
magnetism. 

Silicon carbide has been considered for many years 
as a compound with important potential practical ap- 
plications like power electronics, heterogeneous catalysis 
support, structural and protective components for use 
in future nuclear fusion reactors, etc. This interest has 
been fostered in the last years by the possibility to ob- 
tain SiC in different atomic configurations ranging from 
macro molecules like fullerenes to two-dimensional sheets 
and their wrapped configuration in nanotubes (see for 
instance the work of Melinon et al. pjQj and references 
therein) . 

The free SiC (110) surface has been already theoret- 
ically studied by Sabisch et al. [TT] (for a complete re- 
view of the different SiC surfaces see reference [H]). Un- 
like other surface orientations, the free SiC(llO) surface 
does not present any reconstruction as different theoreti- 
cal works reveal [12] . To the best of our knowledge there 
are no published experimental studies of this surface. 

Here, in this work, we present the results of a first 
principles calculation of different defects at the SiC(llO) 
surface. The calculations throughout this work were per- 



formed within the density functional theory (DFT)[T3j. 
using the generalized gradient approximation (GGA) [T3] 
for the exchange and correlation. The calculations were 
obtained with the SIESTA [T5], [TB] method, which uses a 
basis of numerical atomic orbitals [17] and separable [18] 
norm conserving pseudopotentials |19) with partial core 
corrections [2D] . We have found satisfactory the standard 
double-^ basis with polarization orbitals (DZP) which 
has been used throughout this work. The bulk calcula- 
tion yields a lattice constant of 4.41 A (Si-C bond length 
of 1.91 A) in fair agreement with the experimental value 
of 4.36 A. We find a partially ionic Si-C bond with a 
charge transfer of 0.45 electrons to the Si atom as eval- 
uated using Mulliken population analysis. We obtain a 
band gap of 1.31 eV similar the 1.29 eV calculated by 
Sabisch et al. [H] and to be compared with the experi- 
mental one of 2.417 eV [2"T] . 

To simulate the free SiC(llO) surface we have per- 
formed calculations of different size slabs formed by the 
stacking of (110) planes. In particular, we have consid- 
ered in this work slabs built by 7 and 15 layers. In all 
cases the "back surface was properly saturated with hy- 
drogen atoms. The convergence of the relevant precision 
parameters was carefully checked. The real space inte- 
gration grid had a cut-off of 400 Ryd. Of the order of up 
to 600 k points were used in the two-dimensional Bril- 
louin zone sampling using the Monkhorst-Pack k-points 
sampling. To accelerate the self-consistency convergence, 
a polynomial broadening of the energy levels was per- 
formed using the method of Methfessel and Paxton [55] 
which is very suitable for systems with a large variation 
of the density of states at the vicinity of the Fermi level 
which is the case in our system (see below). Broadening 
like Fermi-Dirac can be inappropriate and give wrong 
results. It is worth to mention that the energy differ- 
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FIG. 2: (Color online) Projected electronic densities of states 
at the surface C (solid black line) and Si (broken blue line) 
atoms in the SiC(llO) surface. The origin of energies is taken 
at the Fermi Energy. A small imaginary contribution to the 
energy has been added for presentation purposes. 



FIG. 1: Calculated band structure an total density of states 
for a 15-layers slab of SiC(llO) with a buckled, unrecon- 
structed clean surface. The surface states band structure is 
apparent at the bulk band gap. The origin of energies is taken 
at the Fermi Energy. 



ences between non-magnetic and magnetic solutions are, 
in general, small, what requires a very high convergence 
in all precision parameters and tolerances. To obtain 
the equilibrium geometry we relaxed the atoms until the 
forces acting on them were smaller than 0.01 eV/A. 

Results of the calculations of a 15 layers slab are re- 
ported in Figures [l] and [2] We, like in reference [TT] . 
obtain a non-metallic surface with two well separated 
surface bands (see Figure [lj ; the occupied narrow one 
is mainly due to the C dangling bond , whereas the un- 
occupied wider one has a Si character as the projected 
densities of states reveals (see Figure[2|. We, like Sabisch 
et al. [11] , obtain a buckled non reconstructed surface; 
the C and Si atoms are displaced normal to the surface 
0.01 A outwards and 0.24 A inwards respectively . The 
surface bonds are shorter (1.81 A) than the bulk ones 
(1.91 A) and the bond angles (~ 122°) at the surface 
are larger than at the bulk (109.47°) making the sur- 
face chains more one-dimensional than the ideal structure 



ones. We allowed in the calculation the possibility of dif- 
ferent spin populations in the atoms. We did not find 
a stable solution with non-zero magnetic moment at the 
surface atoms. It is worth to indicate that the densities 
of states at Figure [2] suggest that introduction of holes 
in the system would shift the Fermi level to the sharp C- 
like peak giving rise eventually to a removal of states at 
the Fermi level by means of structural (Jahn- Teller type) 
and/or magnetic broken symmetries. However, introduc- 
tion of electrons would shift the Fermi level to a much 
wider Si-like band, magnetic solutions being unlikely. 

We next consider different defects at the surface. We 
have considered a 2x2 surface unit cell with one defect in 
it. Nine (110) layers were included in the slab calcula- 
tion. In most cases the results were checked by enlarging 
the supercell size. We first study a hydrogen atom sat- 
urating either a Si or a C surface atom. In the case 
of hydrogen bonded to a Si atom we obtain, after full 
geometry relaxation, a magnetic solution such that the 
surface carbon atoms next to the defect have a magnetic 
moments of 0.705 \ib and 0.058 [Ib respectively (see Ta- 
ble [i}. This magnetic moment is very localized at one of 
the carbon atoms, the other carbon atoms in the unit cell 
have a magnetic moment of the order of 0.02 fis whereas 
the induced magnetic moments at the silicon atoms is 
also very small (see Table [T]) . It should be indicated that 
the calculated total magnetic moment in the slab is 1 fis ■ 
This magnetic solution is the stablest one with an energy 
gain of 0.142 eV per unit cell with respect to the diamag- 
netic one. The spin resolved projected densities of states 
are shown in Figure[3j The splitting of the carbon surface 
state peak due to the local magnetic moment is apparent. 
It is interesting to realize that the magnetic moment sym- 
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FIG. 3: (Color online) Projected spin resolved densities of 
states at the vicinity of an hydrogen atom chemisorbed at 
the SiC(llO) surface. Solid (broken) lines represent spin up 
(down) states. Panels (a) and (b) stand for the densities of 
states projected on the C atom closest to the defect (H bonded 
to the surface Si atom). Panel (c) shows the densities of states 
projected on the Si atom closest to the defect (H bonded to 
the surface C atom) . The arrow at the insets indicate the 
atom where the density of states is projected. The origin of 
energies is taken at the Fermi Energy. A small imaginary 
contribution to the energy has been added for presentation 
purposes. 



(a) H at the SiC-(1 1 0) Surface 
Spin-independent Solution 
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FIG. 4: Projected densities of states near adsorbed hydrogen 
in the spin independent solution, (a) Projected density of 
states at the Si surface atom next to the carbon atom bonded 
to the hydrogen impurity. (b)Projected density of states at 
the C surface atom next to the silicon atom bonded to the 
hydrogen impurity. In both cases the origin of energies is 
taken at the Fermi Energy 



metry of the two C atoms bonded to the defect is broken 
(see Table [I]) and, in addition, the symmetric geometry 
is also broken, including a buckling of 0.22 A of the C 
atoms. It is worth to note that while the magnetic solu- 
tion is accompanied by this Jahn- Teller type distortion, 
a spin independent calculation does not give rise to any 
structural distortion. Our result clearly indicates an in- 
teresting lattice distortion local magnetic coupling. The 
same calculations were performed in the case of H at- 
tached at the surface carbon atom. In this case, as antic- 
ipated, no magnetic stable solution was found obtaining 
a paramagnetic ground state. 

In order to understand the origin of these localized 
moments when H is bonded to Si, and the absence of 
magnetic moments at the Si atoms when H is attached 
to C, we present in Figure [4] the results of the projected 



densities of a spin independent calculation. We first ob- 
serve that in both cases the Fermi level lies at a band 
and therefore the system is metallic. However, the main 
difference between C and Si is that the density of states 
at the Fermi level in the case of H on top of Si is much 
higher than in the case of H on top of C. The C band 
is much narrower than the Si band due to the more lo- 
calized valence orbitals in C than in Si, the correspond- 
ing covalent radii are 1.17 A and 0.77 A for Si and C 
respectively. This local character in the C case is also 
responsible for a larger electron-electron Hubbard like in- 
teraction U and therefore the fulfillment of the condition 
U x N(Ep) > 1, with N(Ep) the substrate density of 
states at the Fermi energy, to develop localized magnetic 
moments in a Stoner-like approach |23| . 

Other possible surface defects were considered in de- 
tail: 
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TABLE I: Charge (Q) and magnetic moment (n) at carbon 
atoms near different defects at the SiC(llO) surface. The 
charges are the Mulliken populations. Q f,i stand for the 
spin up and spin down electronic charge respectively. 

Atom Q t Q I Qtot fi (hb) Si-C (A) 
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C next to extra C 2.071 1.855 3.926 0.216 



C next to extra C 2.071 1.855 3.926 0.216 
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i) A fluorine impurity. F behaves like H, it satu- 
rates the Si dangling bond and induces magnetism in 
the nearby C atoms (see Table [T] and Figure [5jc)). The 
magnetic moments and atomic relaxations induced by F 
are very similar to those induced by H. Like in the case 
of H the total magnetic moment induced in the crystal is 

1 fJL B . 

ii) An oxygen impurity. Oxygen can saturate the sur- 
face dangling bonds but energetically is much more fa- 
vorable to break the Si-C surface bond rather than satu- 
rating the dangling bonds, like in the formation of Si02 
at the Si surface [53]. In this case, no magnetic solution 
was found. The corresponding projected density of states 
is drawn in Figure [5jb). 

hi) A substitutional nitrogen atom. The most stable 
configuration of an N impurity is by substituting a Si 
surface atom and being bonded to three C atoms. The 
results of the calculation of this defect display no mag- 
netic moment. The density of states at one of the surface 
carbon atoms bonded to N is shown in Figure [5|a). 

iv) A silicon vacancy. This configuration favors a mag- 
netic solution at the C atoms (see Table [I]) such that the 
magnetic moment at the underneath carbon atom is in 
antifcrromagnetic configuration with respect to the two 
surface C atoms. The total magnetic moment in the slab 
is zero. The projected densities of states at the carbon 
atoms close to the vacancy are shown in Figure [6] 

v) Substitutional carbon at the surface layer. In this 
case we also obtain a magnetic solution of the carbon 
atoms near the defect. The magnetic configuration of 



FIG. 5: (Color online) Projected densities of states near dif- 
ferent impurities at the SiC(llO) surface. The impurities are 
indicated by a black circle. The arrows indicate the C atoms 
where the density of states is projected. In the case of flu- 
orine (panel (c)) both spin up (continuous black) and spin 
down (broken red) densities of states are shown. In the case 
of oxygen (panel (b)) and nitrogen (panel (a)) there is no dis- 
tinction between spin up and spin down bands. The atomic 
configurations at the different defect geometries are sketched. 
The broken lines in panel (a) indicate the carbon atom under 
the topmost layer 

the carbon atoms is antiferromagnetic like in the case of 
the C-(110) surface [§]■ The total magnetic moment in 
the slab being also zero. 

In conclusion, we have found that different defects at 
the SiC(llO) surface can induce localized magnetic mo- 
ments at the C atoms near the defect with various mag- 
netic configurations. In the case of hydrogen or fluorine 
adsorbed at a Si surface atom magnetism is accompa- 
nied by a local lattice distortion. Magnetism takes place 
mainly at the carbon atoms due to the localized charac- 
ter of its atomic valence charge. We have found that the 
SiC(llO) surface is an excellent candidate for magnetic 
configurations due to the quasi one dimensional behavior 
of the topmost Si-C chains. 

We thank J. M. Soler for helpful discussions about this 
work. Financial support of the Spanish Ministry of Sci- 
ence and Innovation through grants FIS2009-12712 and 
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FIG. 6: (Color online) Projected densities of states near a Si 
vacancy at the SiC(llO) surface. The surface C atoms are 
indicated by circles in the sketches. The broken lines in the 
sketches indicate the carbon atom under the topmost layer. 
The C atoms close the vacancy where the densities of states 
are projected are indicated by an arrow, (a) Densities of states 
at the outermost C atoms, (b) Densities of states at the C 
atom at the underneath layer. Solid black (broken red) lines 
stand for spin up (down) densities of states. 
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